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In view of ongoing experiments to trap ultracold spin-polarized 6 Li, we study various properties of 
an interacting Fermi gas in a harmonic trap taking the discrete nature of the unperturbed harmonic 
00 ■ trap levels into account exactly. As 6 Li has a rather large and negative scattering length, we focus on 

Q\ | the effects of the attractive atom-atom interaction on several thermodynamic properties and on the 

momentum and density distributions. The dependence of the chemical potential, the specific heat 
and the density and momentum distributions on the number of particles in the trap is obtained. We 
also calculate the energy of the gas. Comparison is made with results of a semiclassical calculation 
and with the properties of a non-interacting gas. We find that the effect of the interactions is 
rather large for realistic trap frequencies. Hence, it is important to include these interactions in any 
quantitative predictions relevant for experiments. 

<N 

T^ 1 I. INTRODUCTION 

a: 

Recently there has been a lot of interest in the properties of trapped ultra-cold atoms. This interest has largely been 
sparked by the achievement of Bose-Einstein condensation in the Bosonic systems 87 Rb, 7 Li, and 23 Na Several 
theoretical studies of trapped degenerate Fermi gases have been presented in addition to those on Bose systems. It 
has even been shown that a two-component gas of spin-polarized atomic 6 Li becomes superfluid at experimentally 
obtainable densities and temperatures |Q. This is due to the fact that 6 Li has an anomalously large and negative 
s-wave scattering length a, with a recent measurement giving a ~ — 2160ao where ao is the Bohr radius j^]. s-wave 
c/3 , scattering is forbidden for two fermions in the same spin-state, but, the 6 Li atom has six hyperfine states. This means 
that by trapping 6 Li in two different hyperfine states, one can observe the relatively strong interactions due to s-wave 
scattering between atoms in different states. The trapping can be achieved since the energy of the hyperfine states 
depend on the external magnetic field. The highest three states "prefer" a lower magnetic field and can therefore be 
trapped in a static magnetic trap. In particular, one can trap the two highest states as proposed by Houbiers et al [0. 
These two states will here be labeled | |) and \ [). In such a two-component gas of trapped spin-polarized 6 Li atoms, 
there will only be, to a good approximation, interactions between atoms in different hyperfine states, whereas there 
, will be almost no interaction between atoms in the same state. Indeed, it was for such a two-component system that 
the relatively high transition temperature T c for a BCS-type phase transition was predicted 
£ — , The purpose of this paper is to examine the normal state properties of such a trapped gas of 6 Li atoms with 
components in the two highest hyperfine states. It is clearly necessary to understand these properties in the path 
to achieving the predicted BCS-type transition. The normal state properties of a non-interacting trapped cloud of 
fermions have been treated within the semiclassical Thomas-Fermi approximation and also by taking the discrete 
nature of the trap energy levels in the non-interacting limit into account ||. As the interactions between 6 Li atoms 
in two different hyperfine states are relatively strong, it is important to include the effect of these interactions in any 
realistic treatment of the system. Hence, the present paper extends the analysis of Refs. by including both the 
discrete nature of the trap levels as well as the effects of the interactions (in the mean-field approximation). 
■ The paper is organized as follows. In Sec. || we set up the formalism needed to treat an interacting two-component 
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system of fermions within the mean-field approximation. We then, in Sec. Ill, analyze the influence of the trap 
potential and the interactions on the chemical potential of the gas. We show, that the interactions have two distinct 
effects on the chemical potential: they lower its value below the non-interacting value and smooth out the steplike 
features predicted by Schneider and Wallis Q . The quasiparticle energy spectra and wave functions are considered in 
Sec. |^ and in Sec. |v| we investigate the behaviour of the energy and the heat capacity. We find that at sufficiently 
low temperatures the interactions change the qualitative behaviour of the heat capacity. This is explained in terms 



of the effect the interactions have on the quasiparticle spectrum. In Sec. VI we discuss the density and momentum 



distribution and their deviati ons f rom the non-interacting case and from the Thomas-Fermi predictions. Finally, we 



summarize the results in Sec. VII 
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II. FORMALISM 



We consider a dilute gas of interacting 6 Li atoms in two hyperfine states trapped in an external potential Uq(t). As 
the gas is dilute the interactions mainly happen through two-body collisions. Furthermore, since the s-wave scattering 
length is much larger than the p-wave scattering length we can neglect any interaction between fermions in the same 
hyperfine state. The gas is then described by the Hamiltonian: 



II=J2[ d3r V4(r)[^V 2 + imwV - /# ff (r) - gj d 3 r V>f(r)^(r)^(r)V> T (r), (1) 
(j 

where m is the mass of the particles, the attractive interparticle potential has been approximated by a contact 
potential V(r' — r) ~ ~gS(r' — r), g > 0. The field operators ip a (r) obey the usual fermionic anticommutation rules 
and describe the annihilation of a fermion at position r in the hyperfine state \a). The trapping potential is for 
simplicity assumed to be well described by an isotropic harmonic oscillator Uq(t) — ^mu> 2 r 2 . The trapping frequency 
u> is taken to be the same for each hyperfine state. We have assumed that number of particles N in each state is 
the same such that we only have one chemical potential /i. As the critical temperature for a BCS type transition is 
maximum when the number of particles in the two hyperfine levels is equal , we expect this configuration to have 
the most experimental relevance. The non-interacting case is achieved by setting g — 0; this limit has been treated by 
Butts and Rokhsar within the Thomas-Fermi approximation and by Schneider and Wallis ^] taking the quantizing 
effect of the trap potential into account. In this paper, we are interested in the effect of the interactions on the normal 
state properties of the gas. We can, therefore, ignore any pairing correlations leading to a BCS type transition and 
use the following mean-field Hamiltonian: 

H M ean = £ / ^^Mh^V 2 + + U (r) - /#«r(r). (2) 

<7 

Here the self-consistent field U(r) = — <?(^>l .(r)'i/>_ cr (r)) is the standard Hartree-Fock result for a contact interaction. 
To include all two-body scattering processes on the mean-field level, one can put g = 4Tr\a\h 2 /m, with a being 
the s-wave scattering length of collisions between the fermions in the two hyperfine states ||. As the trapping 
potential is assumed to be isotropic, the self-consistent solution of the mean-field Hamiltonian with the lowest energy 
will be spherically symmetric. The Hamiltonian can then readily be diagonalized by writing the field operator as 
ip a (r) = J2vim u im( r ) a ahn ■ The operator a^ lm describes the annihilation of a quasiparticle in the hyperfine state a 
with total angular momentum [1(1 + l)^ 2 } 1 ^ 2 , a component along an arbitrary z-direction of rah and a wave function 
u\ m (r). We write the quasiparticle wave function in the form 

ii v (r\ 

uUr) = —Yi m (e,4>), (3) 
r 

where Yi m (9, <fi) are the usual orbital angular momentum eigenfunctions. HMean is then diagonalized by solving 

= + ^Z^ 1 + \™Jr 2 + U(r) „K(r) (4) 



for each angular momentum I. Here the quasiparticle energy is independent of m due to the spherical symmetry. 
The self-consistent potential is determined by 

U(r) = - ff ^|< m (r)| 2 /(^) = -flE *? 2 21 tj ^l) ( 5 ) 

vim ul 

where f(x) = l/[exp(x/ksT) + 1] is the Fermi function and the addition theorem for spherical harmonics has been 
used. Numerically, we use a cutoff of the order E" ^ 2/i for the sum given in Eq. (g). This is more than sufficient since 
we in this paper are considering low temperatures (ksT < Tiuj) and we know that levels with energies higher than a 
few UbT give a vanishing contribution to the density. Thus, within the mean-field approximation, the problem of an 
interacting gas of fermions in two hyperfine states in an isotropic harmonic trap is equivalent to solving Eq. (Q)-(|^) 
self-consistently. Once a solution is found, one can easily calculate various observables for the gas. The total energy 
E for the particles in both hyperfine states is given by 
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E=(H) + 2[iN 

= 2£(2l + l)/(£n J d 3 r< (r)(^V 2 + im W VK (r)-i / d\U{rf. (6) 
^/ 

Here we have added 2N to (H) as it is the total number of particles in the trap. In the following sections we will 
solve the above equations for various sets of parameters and calculate several observables from the solutions. 

As the negative scattering length introduces an attractive interatomic potential, the system can collapse to a fluid or 
solid state when the density of particles becomes too large. This effect has been examined by Houbiers et al. || within 
the semiclassical Thomas-Fermi approximation. For a trap with an equal number of particles in each hyperfine state 
they find, that the spinoidal point is given by N^ 6 \a\/l ~ 0.66 with I = y/KJ cum being the trap length. Numerically, 
this transition is seen from the fact that there is no self-consistent solution to Eq. when the chemical potential 

is too high for a given coupling strength. The density of particles increases for each iteration without bound indicating 
that the system is collapsing into a new dense phase. We find that this problem arises in the region of parameters 
where A^ 1 / 6 1 a | /Z ~ O(l) in qualitative agreement with Houbiers et al. However, due to the computation load we have 
not been able to verify in detail the prediction N x / % \a\/l ~ 0.66 for the spinoidal transition line. As we do not have 
an appropriate theory for such a phase we will in this paper work with parameters (N, g) such that we are well below 
this spinoidal phase-transition line. 

We shall compare some of our results obtained from the solution of Eq. (Q)-(|5|) with approximate results based on 
the Thomas- Fermi approximation. This approximation essentially treats the trap potential as being locally constant. 
The Hartree-Fock equations are then trivially solved for each r by plane waves. The quasiparticle energies are given 
by 

h 2 k 2 

E k (r) = — - gp(r) - fi{r). (7) 

Here the local chemical potential is fi(r) = fj, — ^mu) 2 r 2 and the density is given by 

f H 3 k 

Eq. (0)-(|j) has to be solved self-consistently at each point r. The total number of particles in a single hyperfine state 
in the system is then N — J d 3 r p(r). 



III. THE CHEMICAL POTENTIAL 



In this section we will determine the chemical potential /i as a function of the number of particles in the trap. 
Schneider and Wallis M have done an extensive analysis of f/,(N) in the case of a non-interacting gas. They found 
some remarkable steplike features in /i(N) as compared to the Thomas-Fermi approximation; these steps were due 
to the shell structure of the energy spectrum of an isotropic harmonic trap. They did, however, predict that for a 
real Fermi gas, the interactions would tend to smooth out the steps. Using the formalism outlined above, we are now 
able to examine in detail how the interactions affect these steplike features. In Fig. (1) we have plotted the chemical 
potential \i as a function of the number of particles N in a single hyperfine state for various coupling strengths g for 
a very low temperature (ksT <C hcu). These curves were obtained by solving Eq. (^)-(^|) for varying u, and g. The 
dashed lines correspond to the Thomas-Fermi approximation obtained from the solution of Eq. In Fig. (1) 

(a) we have compared fi(N) for the rather large coupling strengths g' = 1 and g' = 2 with the g = case when there 
are relatively few particles in the trap in order to highlight the stcplike features. We have defined g' = g / (Tiul 3 ) . As 
can be seen, that the interaction has two effects. First of all, it lowers the value of the chemical potential for a given 
number of particles N as compared to the g = case. This is as expected, since the mutual attraction between the 
particles in the two hyperfine states lowers the energy of the gas. Secondly, we note that the step like features of fJ,(N) 
are smoothed out by the interaction. For g' — 1 the step like features survive up to N 2000 and for g' = 2 they 
survive for ^£500. Fig. (1) (b) show fi(N) for g' = 0.4 and g' = 0.2. Putting \g\ = 4Tr\a\H 2 /m, with a = -2160a 
being the s-wave scattering length for 6 Li, g' — 0.4 corresponds to a trap frequency of v — uj/2tt — 144Hz and g 1 = 0.2 
to a trap frequency of v = u>/2ir — 288Hz. These values are close to the ones used in the Bosc-Einstein condensation 
experiments at JILA and at MIT pi] respectively. We again see the same behaviour as in Fig. (1) (a). The step 
like features break down for N w 4 x lCr for g' — 0.4. For g' = 0.2 they break down at N w 10 5 which is not shown 
in Fig. (1). 
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IV. THE ENERGY SPECTRUM 



From Fig (1) we see, that the step like features survive even though the value of fi(N) is several Twj lower than 
the g = prediction. This is somewhat surprising as the steps are associated with the gaps of size Tilo in the 
g = spectrum. One might expect these gaps to disappear once the levels are lowered by more than hu> due to the 
interactions. To examine this effect, we have in Fig. (2) plotted the quasiparticle spectrum for g' = 2 for a very low 
temperature for each angular momentum I. In Fig. (2) (a) the chemical potential is fi/hu> = 8 whereas fi/tkv = 13 in 
Fig. (2) (b). The x-axis denotes the angular momentum I of the quasiparticle states and the quasiparticle energies 
are marked by the symbol x . For comparison the symbols o denote the non-interacting quasiparticle energies. The 
chemical potential is indicated by a thin dotted line. In the non-interacting case (g = 0), the energy spectrum is given 

by E v = (u+ 3/2)hoj with v = 0, 1, 2 The degeneracy of each level is D v = (y + 1)(V + 2)/2 corresponding to the 

angular momentum I being I = 0, 2, . . . v for v even and I = 1, 3, ... v for v odd. This degeneracy in I gives rise to the 
shell structure of the spectrum which is indicated in Fig. (2) by horizontal dashed lines. In the interacting case, the 
energy E" in general depends on both v and I. From Fig. (2) (a) (/j, = 8ft,u>) we see, that the quasiparticle energies in 
the interacting case are lowered several Tioj as compared to the normal state energies; e.g. the lowest energy is —\.lhw 
for g' = 2 as compared to 1.5hu> for g' — 0. However, from the solid lines in Fig. (2) (a) we see, that the energies 
still depend only weakly on the angular momentum I. This means that the Hartree potential lowers the energy of the 
I = v state by almost the same amount as it lowers the / = [y even) or I = 1 (y odd) state. This is not a priori 
obvious as the wave functions for those two states are completely different. This is shown in Fig. (3) where the solid 
lines denote the quasiparticle wave functions u\=® (r) (a) and ^^(r) (b) for fi — 8hu and g' = 2. These states have 
the energies = 8.2h,w and E[Jr% = 8Ahcu respectively. The I = wave function is spread out from r — to 

r w r c i ass = y // 2E n /muj 2 = Al whereas the / = v wave function is peaked around r rs 2.81. Hence one might expect 
that the Hartree potential, which is also plotted as a solid line in Fig. (3) (c), would affect the quasiparticle states 
completely differently. But in fact, the lowering of the two energies as compared to the g = case is approximately 
the same. This explains the fact depicted in Fig. (1) that even though the chemical potential is several Tiui below 
the g = prediction, fi(N) still exhibits the step like features. The reason is, that even though each quasiparticle 
energy is lowered several huj due to the interaction, there is still an approximate Z-degeneracy: there are still bands 
separated in energy by « hu in the energy spectrum as a function of I. These bands (shell structure) give, as in the 
g = case, rise to the steplike structure of /J,(N). For comparison, we have also in Fig. (3) plotted the non- interacting 
wave functions as dashed lines. As expected, we see that the effect of the interactions is to compress the quasiparticle 
states closer to the center of the trap where the Hartree field is large. The cloud of particles is compressed as will 
also be seen from the density distributions plotted in Sec. (|Vlj). We can qualitatively understand the slight increase 
in energy with increasing I as seen from the solid lines in Fig. (2) (a). The low angular momentum functions have a 
larger amplitude in the center of the trap and are thus more affected by the Hartree field than the high I states. 

Contrary to the /i = 8Huj case, we see from the solid lines in Fig. (2) (b), that when /i — lZfim the approximate 
independence of the quasiparticle energies on I no longer holds. The Hartree field is now so strong that it has washed 
out the shell structure of the spectrum and it is qualitatively different from the g = case. The energies now increase 
significantly with increasing I. This explains the fact from Fig. (1) (a), that for [i ~ 13 and g' = 2 the steplike features 
have disappeared. 



For a gas of particles in a constant confining potential, the most useful definition of the heat capacity is [[L2| is 
C N = ^dE/dTlN with E(T, N) given by Eq. (||). As pointed out by Schneider and Wallis ||, the shell structure of 
the harmonic trap spectrum has a drastic consequence for the low temperature (fc^T <C tuv) heat capacity; the gaps 
in the energy spectrum for g = mean that the heat capacity is exponentially suppressed for low temperatures. In 
the non-interacting case the total energy for particles in both hyperfine states is given by 



When the number of particles is such that a finite number of energy levels up to and including the level E UF are 
completely filled for T — (i.e. N = 1, 4, 10, 20, . . . for vf = 0, 1, 2, 3 . . .), the low temperature chemical potential for 
a constant number of particles is given by 



V. THE ENERGY AND THE HEAT CAPACITY 




(9) 




(10) 
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with k„ f = D VF+ i/ D VF > 1. The low temperature heat capacity is then easily found to be 

C [T\ 1 (hi A^p—Pfiu/z 

_ [,„ + 5/2) D „ p+1 _ - ( W + (") 

where /3 = 1/kgT. This is, as one would expect, suppressed by a factor exp(— f3fkd/2) as compared to the usual low 
temperature Thomas-Fermi result: 



2Afc s /^(GAT) 1 / 3 ' 



(12) 



It is interesting to examine how the interactions change this result. In Fig. (4) we have plotted the energy per 
particle E'/2N of the interacting gas as a function of temperature T' = T/fico, where E' = E/hu and the energy is 
given by Eq. (^|). We have used the parameters g' — 0.4, N = 5456 (a) and g' = 0.4, N = 43680 (b). For comparison 
we have also plotted the g — results as dashed lines in Fig. (4) (c)-(d). As can be seen the interactions for both set 
of parameters have lowered the total energy of the gas considerably. This is as expected as we from Fig. (1) see, that 
the chemical potential in both cases is significantly lower than the g = result. Fig. (5) depicts the corresponding 
heat capacity as a function of the temperature for N — 5456 (a) and N — 43680 (b) (solid lines). The dashed lines are 
the g = results and the dash-dotted line is the Thomas- Fermi result for g' = as given by Eq. (|l2|). For N = 5456 
(Fig. (5) (a)) we see, that the heat capacity for the interating system is still exponentially suppressed for ksT <C Tuo. 
For very low T the heat capacity is practically zero as for the g = case and the interactions have not changed the 
qualitative behaviour of the heat capacity. This is a direct result of the fact that the interactions have not changed 
the shell structure of the energy spectrum. This can be seen from Fig. (1) (b)) where the step like features of p(N) 
are still prevailing for g' = 0.4 and N = 5456. 

The break-down of this effect is shown in Fig. (5) (b). We see that for N = 43680 the interactions have now 
increased heat capacity substantially over the non-interacting result for low temperatures. It has the same order of 
magnitude as the Thomas-Fermi prediction and it is qualitatively different from the non-interacting case. This is due 
to the fact that for this set of parameters the interactions have washed out the shell structure in the energy spectrum 
(the i-degeneracy) as can be seen from Fig. (1), where the steplike features are smoothed out for N = 43680 and 
g' = 0.4. There are no gaps of hui in the spectrum and hence no exponential damping factor in the heat capacity. 

It should be noted that for very low temperatures the possibility of a transition to a superfluid state arises [Q. 
From the well known result of weak coupling superconductors, we expect that this transition will give rise to a kink 
in the E(T) curve as the energy per particle is lowered as compared to the normal state p~3| ] ■ From this kink, we 
obtain a discontinuity in Cn(T) at the transition temperature T c . We will not in this paper consider this effect as 
we are concentrating on the normal state behaviour. Work is under progress to examine for which temperatures and 
densities this transition takes place and how it can be determined experimentally. 

We conclude, that the interactions in general lower the total energy as one would expect. The effect of the interaction 
on the heat capacity depends on the strength of the interactions and on the number of particles in the system. When 
few particles are in the system or when the interaction is so weak that the shell structure of the energy spectrum is 
intact, the heat capacity is still exponentially suppressed for low temperatures. The heat capacity can be suppressed 
although the total energy of the system is significantly lower than for the non-interacting case. This is, of course, 
a consequence of the fact that the lowering of the quasiparticle ener gies and the washing out of the shell structure 



happen at different set of parameters g and N as explained in Sec. Ill -IV. Once the interactions are strong enough to 
wash out the shell structure around the chemical potential, the behaviour of the heat capacity changes qualitatively. 
It is no longer suppressed and has the same order of magnitude as the Thomas-Fermi prediction. The transition 
between those two limiting behaviours is smooth as a function of TV or g. 



VI. DENSITY DISTRIBUTIONS 



In this section we will present results for the density and momentum distributions. In an isotropic trap, these 
distributions will be spherically symmetric when the gas is in the ground state. The density distribution p(r) is 
calculated from p(r) = —U(r)/g with U(r) given by Eq. (||). In Fig. (6), p{r/l) is displayed for N — 120,165 and 
g' = 2. For comparison we have also plotted the distributions for g' = 0. The dashed curves are the Thomas- Fermi 
results obtained from the solution of Eq. (Q)-(||) with a given N. The temperature is taken to be zero. As can be 
seen, p{r/l) is substantially changed due to the interactions. The cloud of particles is compressed as compared to the 
non-interacting result due to the attractive forces. As a high density of particles increases the critical temperature for 
a BCS-type transition, this effect favours the formation of the superfluid state M. Furthermore, we observe a central 
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minimum for TV = 120 and a central peak for TV = 165 for both g' = 2 and g' = as compared to the Thomas-Fermi 
predictions. For the g' = 0, this is because N — 165 corresponds to a filled shell with vp = 8, N = 120 corresponds 
to a filled shell with vp = 7 and the fact that shells with v = odd do not contribute to p(0) as they have odd angular 
momentum We see that the interacting system exhibit the same qualitative dependence of p(0) for this set of 
parameters, even though the actual densities are substationally different from their g' = counterparts. This is due 
to the fact, that for this set of parameters, the shell structure of the quasiparticle spectrum is still intact although the 
actual energies and wave functions are changed substantially. This can be seen from Fig. (1) (a) where the step like 
features still prevail for this set of parameters. Thus, even for g' = 2, N — 165 still corresponds a highest shell of even 
angular momentum states totally filled whereas N = 120 corresponds a filled highest shell of odd angular momentum 
states. For a larger number of particles, where the shell structure of the quasiparticle spectrum has been washed out, 
it turns out that this behaviour of p(0) has disappeared as expected. Also, a non-vanishing temperature will tend to 
wash out the predicted dip/top behaviour of the density as the transition from occupied to unoccupied shells as a 
function of energy becomes less abrupt. 

The momentum distribution of the particles is not completely straightforward to measure. A simple free expansion 
experiment where one switches off the trap potential non-adiabatically does not strictly measure this distribution. 
The reason is, that the gas does not expand freely; collisions between particles in different hyperfine states will alter 
the momentum distribution. One would need a time-dependent formalism |]l4| in order to treat such an expansion 
rigorously. However, we will here assume that the momentum distribution in the trapped state gives a good indication 
of the distribution measured in such a free expansion experiment. The spherically symmetric momentum distribution 
(c\ck), where cl creates a particle in a plane wave state exp(ikz) along an arbitrary z-direction, is calculated from 

{clc k )=J2\(kKo)\ 2 f(En- (13) 

Here we have utilized the fact, that a plane wave along the z— direction only contains m — spherical harmonics. 
Using the well known expansion of a plane wave exp(ikz) in spherical harmonics |l5| we obtain 

<AKo> = / d\e-^u^{v) 

/ ■jr \ 1/2 f°° 

= H)V(2^ + l)47r(-) J J l+1/2 (kr)uUr)V^dr (14) 

where Ji+\/2{x) is the ordinary Bessel function. In Fig. (7), we have plotted the momentum distribution for N — 43680 
and g' = 0.4. The momentum is measured in units of I . The dashed curve is the non-interacting result. For this 
high number of particles, the shell structure of the quasiparticle spectrum is washed out and the distributions are, 
apart from a small shell around the edge of the cloud, almost identical to the Thomas-Fermi prediction (not plotted) . 
We see that the interactions have spread out the distribution considerably as compared to the non-interacting case. As 
the Hartree field lowers the quasiparticle energies more levels become occupied in the center of the trap. The higher 
momentum states thus become populated leading to a spreading out of the momentum distribution as compared to 
the g' = case. 

We conclude that the interactions in general alter both the momentum distribution and the density distribution 
substantially. The density distribution is compressed and the momentum distribution spread out as compared to the 
non-interacting results. Furthermore, the central minima and maxima of p(r) as a function of N can still be observed, 
when the interactions have not yet washed out the shell structure of the quasiparticle spectrum. 

VII. CONCLUSION 

In this paper we have considered a trapped spin-polarized gas of interacting fermions. We find that the interactions 
have two distinct effects on the quasiparticle spectrum. It lowers the quasiparticle energies and thus the total energy 
of the gas. Also, above a certain number of particles in the trap, it washes out the shell like structure of the spectrum 
associated with a harmonic trap. These two effects are independent, in the sense that the energy of the gas can be 
lowered considerably as compared to the non-interacting case, but the shell structure of the spectrum is left relatively 
intact. One can still, therefore, observe effects such as step like features in the chemical potential, exponential 
damping of the low temperature heat capacity and maxima and minima in p(0) associated with a non-interacting 
gas. Whether these effects will be observable depend on whether the condition ksT <C fiio is experimentally feasible. 
The interactions are also found to compress the atom cloud and spread out the momentum distribution considerably. 
This effect should be readily observable; it is important to include in any realistic calculation of the properties of 
spin-polarized 6 Li in a trap. 
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Figure Captions 



Fig. 1: /u(TV CT ) in units of Tito for g' = 1, g' = 2 (a) and </' = 0.2, g' = 0.4 (b). The dashed lines are the Thomas-Fermi 
approximation. 

Fig. 2: The quasiparticle spectrum in units of Tito for g' = 2 , ^i/Tkj = 8 (a) and g' = 2, £i/?iu> = 13 (b). 

Fig. 3: The radial quasiparticle wave functions for v = 8, I = (a) and v = 8, I = 8 (b) as a function of r/l. The 
solid lines arc for g 1 = 2 and the dashed lines are for g — 0. The Hartree field in units of Tito is plotted in (c). 

Fig. 4: The energy in units of Tiui for g 1 = 0.4 (solid lines) and TV = 5456 (a) and N = 43680 (b). The dashed lines in 
are for the non- interacting case for N = 5456 (c) and TV = 43680 (d). 

Fig. 5: The heat capacity in units of k B for g' = 0.4 and TV = 5456 (a) and TV = 43680 (b). 

Fig. 6: The density distribution p(r/Z)for TV = 120, 165, g' — 0, 2 (solid lines). The dashed lines are the Thomas-Fermi 
results. 

Fig. 7: The momentum distribution (c\c k ) for TV = 43680, g' = 0.4 (solid) and TV = 43680, g' = 0.0 (dashed) 
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